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Abstract 

We study a hermitian (n + l)-matrix model with plaquette interaction, Y!2=i MAiMAi. 
By means of a conformal transformation we rewrite the model as an 0(n) model on 
a random lattice with a non polynomial potential. This allows us to solve the model 
exactly. We investigate the critical properties of the plaquette model and find that for 
n g] — 2,2] the model belongs to the same universality class as the 0(n) model on a 
random lattice. 
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1 Introduction 



Despite the fact that great progress has been made in solving matrix models in re- 
cent years many interesting models remain unsolved. One important class of models 
for which an exact solution is still lacking is models with "plaquette type" interac- 
tions. Lattice gauge theories like the Weingarten model and the Kazakov-Migdal 
model are typical examples of such models but recently also plaquette type models 
without gauge degrees of freedom have attracted attention, namely as generating func- 
tional for Meander numbers M. In the present paper we will consider the following 
model 



N 2 F 



JdAffldAiexpl-NtT (v(M) + ^A 2 i -^MA i MAM (1.1) 

i=l (. \ i=l i=l I ) 



where all the matrices are hermitian and V(M) is an arbitrary polynomial potential. 
When n = 1 this model shows a large degree of similarity with the 2-dimensional 
reduced Weingarten model |l], U which is given by 



Z = J II d^dA^exp | -iVtr 
n=i { 



£ 4 A, - g £ (A,A v A\Al + h.c 

[1=1 fl,U=l 



1.2) 



However, the two models are not equivalent. A model equivalent to ( |1 . 1| ) for n 
involving complex matrices is 

Z = J dM 1 dM 2 dA^dAexp j-JVtr V{M X ) + ^M 2 2 + A ] A - M X .\'M,.\ 



where M\ and M 2 are hermitian and A is complex. Our model (|1 . 1|) also shows some 
similarity with matrix models generating Meander numbers M. Its interaction is of the 
type needed for such models. However, our model is too simple to provide a generating 
functional for Meander numbers. For that purpose one must be able to work also with 
an arbitrary number of M-matrices. Let us finish by mentioning that our solution 
of the model (1.1) gives the solution to a simple three-matrix problem, namely the 
following 

Z = J dAdBdCexp j-iVtr (v(A) + ^B 2 + ^C 2 - gABC^ . (1.3) 

The partition function (|1.3| ) can be brought on the form fll.lj) (with n = 1) by integrat- 
ing out one of the three matrices. In reference 0] the model ( |1.3| ) with V(A) = \A 2 
was studied numerically. 

The paper is organized as follows. In section ^| we derive the saddle point equation 
of the model ( |1.1| ) and argue that it has the same structure as that of the 0{n) model 



2 



on a random lattice [EJ. Then in section ^| we explicitly transform the model into an 
0(n) model with a somewhat unconventional potential. Exploiting the already known 
exact solution of the 0(n) model on a random lattice |7j], we hereafter in section [| 
write down the solution of the present model. In section ^| we specialize to a quadratic 
potential and perform a detailed analysis of this case. In particular we investigate the 
critical properties of the model and find that for n e] — 2, 2] the model ( |1.1|) belongs to 
the same universality class as the ordinary 0(n) model on a random lattice. Section |] 
contains our conclusion and outlook. Finally in an appendix we comment on the 
Virasoro algebra structure carried by our model. 



2 The saddle point equation 

Let us carry out the gaussian integration over the A- matrices in (p..l|) . This gives 

Z = J dMexp {-Ntr V(M)} det (l <g> I - M T ® m)~™ /2 . (2.1) 

where M T is the transpose of M. Next, let us diagonalize the M- matrices and integrate 
out the angular degrees of freedom. This leaves us with the following integral over the 
eigenvalues, {A;}, of the matrix M 



Z oc 



/ J] dXi exp \-Nj2V (A,) 1 I] (Aj - A fc ) 2 11(1 - \j\ k )~ n/2 . (2.2) 

i { j ) j<k j,k 



In the limit iV — > oo the eigenvalue configuration is determined by the saddle point of 
the integral above ||. The corresponding saddle point equation reads 

NV (A,) = 2 £ - J— + n £ -\- (2.3) 



or 



Following fH we now introduce an eigenvalue density p(A) = J2% $ (A — Aj) which in 
the limit N — > oo becomes a continuous function. As is clear from equation ( j2.2|) the 
model becomes singular if one of the eigenvalues approaches +1 or —1. We shall hence 
solve the model with the requirement that the support of the eigenvalue distribution 
does not include these points. To be precise we will assume that the eigenvalues are 
confined to one interval, [a, 0\ , — 1 < a < (3 < 1 and that the corresponding eigenvalue 
distribution is normalized to one. Using the hence obtained solution we will afterwards 
investigate what happens when, say, a approaches 1. Again following M we introduce 
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the one-loop correlator Q(p) by 

fi(p) = t 'dfi-^-. (2.5) 

J a p — [A 

In terms of the one-loop correlator the saddle point equation ( |2.4| ) can be written as 

n 11 
V'(p) + - = n(p + iO) + {l(p-iO)+n — fi(-), pG[a,/3] (2.6) 
p p A p 

or with V"'(p) = pV'(p) + n and f2(p) = pfi(p) 

H(p + zO) + H(p-zO) + nTt(-) = V'(p), pe[a,/3\. (2.7) 

This equation, in analogy with the saddle point equation of the 0(n) model on a 
random lattice, involves two cuts. First there is the cut of the function Q(p). This 
cut is the physical cut, i.e. the support of the eigenvalue distribution corresponding 
to the matrix M. In addition to the physical cut another cut turns up in the saddle 
point equation, namely the cut of the function fi(-). The singular behaviour referred 
to above corresponds to the situation where the two cuts merge. 



3 Transformation to 0{n) model 

In order to fully exploit the similarity of our model with the 0(n) model on a random 
lattice we will explicitly bring it on the 0(n) model form. For that purpose, let us 
perform the following redefinitions of our matrix fields 

M^i^, ^-+i(l + X) 1/2 ^(l + X) 1/2 (3.1) 

i.e., 

dM -> [det(l + X)Y 2N dX, dAi -> [det(l + X)] N dS*. (3.2) 
Inserting these expressions in our original partition function we get 

ZocJdX JJ dSi exp | - JV tr (v (j^) + (2 - n) log(l +X)+X y £S^jY 

(3.3) 

This model is nothing but the 0(n) model on a random lattice with the somewhat 
unconventional potential 

U(p) = V (j^j + ( 2 " n ) +P) = U °(P) + ( 2 - n ) + P)- ( 3 - 4 ) 
The saddle point equation for this model reads 

W(p + iO) + W(p - iO) + nW(-p) = U'(p), (3.5) 
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where W(p) is the one-loop correlator of the X-field. The physical cut now extends 
from a = (1 + ot)/(l — a) to b = (1 + j3)/(l — (3) and the unphysical cut is the mirror 
image with respect to zero of the physical cut. We note that the point p = 1 lies on 
the physical cut and the point p = —1 on the unphysical cut. As before we expect 
some kind of singularity to occur when the physical and the unphysical cut merge and 
we will have to assume that a > 0. Since the potential includes a logarithmic term we 
might also expect a Penner like singularity to appear, i.e. a singularity corresponding 
to the physical and the unphysical cut degenerating to respectively the point p = +1 
and the point p — — 1 [|l(J . 



4 The solution 

In [|7| an exact contour integral representation of the 1-loop correlator of the 0(n) model 
on a random lattice was written down. The derivation of the exact formula was based 
on the assumption of the potential of the model being polynomial. However, it is easy 
to convince oneself that the formulas remain valid (when written in the appropriate 
way) as long as the potential or rather its derivative does not have any singularities 
which intervene with the physical cut of the one-loop correlator. For our model the 
only singular point of U'(p) is p = —1 (cf. to equation ( p.4|) ) and, as argued earlier, 
this point always lies on the unphysical cut. Hence we can take over the solution of the 
0(n) model on a random lattice from |7). Let us remind the reader of the structure of 
the solution. First, it is convenient to decompose the 1-loop correlator into a regular 
part, W r (p), having no cut, and a singular part, W s {p) 

W{p) = W r (p) - W s (p). (4.1) 

It follows from equation (|3.5| ) that W r (p) is given by 

4 — n z 

while W s (p) is a solution to the homogeneous saddle point equation. As shown in ref- 
erence J7| any solution of the homogeneous saddle point equation can be parametrized 
in terms of two auxiliary functions, G(p) and G(p). More precisely, any such solution, 
S(p), can be written as 

S(p) = A(p)G(p) + B(p) P G(p) (4.3) 

where A(p) and B(p) are regular but not necessarily entire functions. The function 
G(p) is defined by the following three requirements 
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1. G(p) is a solution of the homogeneous saddle point equation corresponding to 
n = 2 cos(z/7r), i.e., 



G(p + iO) + G(p - zO) + nG(-p) = 0. 

2. G(p) is analytic in the complex plane except for the cut [a, b] and behaves as 
(j> — a)~ l l 2 and (p — b)^ 1 ^ 2 in the vicinity of a and b. 

3. G(p) has the following asymptotic behaviour 



G{p) 



2 cos(z/7r/2) p' 



P 



oo. 



These three requirements are enough to determine G(p) uniquely and a completely 
explicit expression for G(p) in terms of theta-functions can be written down [Q. We 
shall not need the detailed form of G(p) for the following but let us mention that, as is 
obvious from the definition, G(p) does not contain any explicit reference to the matrix 
model coupling constants. Furthermore the dependence of G(p) on n appears explicitly 
only via a parameter e given by 



e = a sn (z(l — v)K\ k) , 



k-~ 



(4.4) 



The function G(p) is defined in a way analogously to G(p). Only v is replaced by 1 — v. 
Hence G(p) is a solution of the homogeneous saddle point equation with n replaced by 
—n. Now, if G(p) is a solution of the homogeneous saddle point equation corresponding 
to — n then obviously pG(p) is a solution to the original saddle point equation. This 
explains the appearance of the combination pG(p) in relation (|4.3|) . In a compact form 
the full one- loop correlator, W(p), can be written as (cf. to M) 



W{p) 



A-n 2 



ip G(p) 



duj U'(cu) 



-G( P ) I 



lc 2ni p 2 — lu 2 



p 2 - e 2 ) iG(uj) + ^-uG{u) 



c 2m p 2 — cl^ 



(p 2 -i 2 )uG(co)+p 2 ^G(cu 



(4.5) 



where the contour C encircles the physical cut [a, b] but not the points u = ±p and 
where y^e is defined by 



yfl = J(e 2 - a 2 ) (e 2 - b 2 ) = -ab cn (i(l - v)K') dn (i(l - v)K') 



(4.6) 



Here and in the following we will use the convention that tilded quantities appear from 
untilded ones by the replacement v — *■ 1 — v. If one wants to evaluate the contour 
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integral ( j4.5[) in a specific case, the most convenient line of action is to deform the 
contour into several different contours encircling respectively the points uj = ±p and 
the various singularities of U'(p). The contribution from the poles u = ±p then gives 
rise to the regular part of W(p) while the contribution from singularities of U'(p) gives 
the singular part of W{p). The expression ( ^4.5| ) must be supplemented by a set of 
boundary equations which determine the endpoints of the physical cut, a and h. These 
equations read 

£^VV)GM = 0, (4.7) 

a = jg; («) 

and ensure the correct asymptotic behaviour of the one-loop correlator, namely W(p) ~ 
1/p as p — > oo. In @ it was shown that for the ordinary O(n) model on a random lattice 
the higher genera contributions to the correlators and the free energy simplify consid- 
erably if one expresses the p-dependence via a set of basis functions ^G^(p), Gj®(p)\ 
and the dependence on the coupling constants via a set of moment variables {Mk, J^}. 
Needless to say that a similar simplification can be obtained in the present case. 



5 The quadratic potential 

For simplicity, let us now restrict ourselves to the case where the potential V(M) 
in (|1.1| ) is given by 

V(M) = ^M 2 . (5.1) 

The analysis of the general case can be done along the same lines. For U(p) in equa- 
tion ( |3.4|) we then obviously have 

U (P) = ^(jjfj +(2-n)log(l+p). (5.2) 

5.1 The boundary equations 

Inserting (|5.2|) into the boundary equation ( |4.7|) we get 

(2 - n)G(-l) - - 

c 

Here the first term comes from the pole at w = — 1 in the logarithmic term and the 
second from the pole at u = — 1 in U' Q {p). There is no contribution from infinity. Next 
inserting the expression for U'(p) into (|4.8|) we get 

= 0, (5.4) 

p=-i 



cP_ d_ 

dp 2 dp 



G(p) 



(5.3) 



(2 - n)G(-l) + ~G(p) 
c op 



15=-] 



2 

c 



d 2 
dp' 



+ 1 



G(p) 
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where the first term comes from the pole at to = —1 in the logarithmic term and the 
two next from the pole at uj = — 1 in U' Q {p). In this case we do have a contribution 
from infinity but it cancels with the constant on the right hand side of the original 
equation. To proceed we need to know -j^G(p) and J^-G(p). These can be found 
by exploiting the fact that any solution to the homogeneous saddle point equation 
corresponding to n = 2 cos(^7r) has a parametrization of the form ( f4.3|) (and similarly 
for the saddle point equation corresponding to n = 2cos((l — v)tt)). The exact form 
of the parametrization is determined by the requirements on the analyticity properties 
and the asymptotic behaviour of the functions in question (cf. to 0). The result reads 



dp ^ (p 2 — a 2 )(p 2 — b 2 ) 



p\e 2 -p 2 - a— j G(p) + i (ap 2 + e 2 d) G(p) | , 



(5.5) 



d 2 1 



X 



2/ + I a 2 + b 2 - 5 ( e 2 - I - aa \ p 4 



+ ((a 2 + b 2 ) ^2e 2 -a^-« 2 j 



4a 2 6 2 ] p 2 + a 2 b 2 ( e 2 - ad - a^- 



G{p) 



+ 



4ap 4 + (6e 2 a - a(a 2 + b 2 )) p 2 - 3e 2 d (a 2 + b 2 ) - 2aa 2 b 2 ] (-i)pG(p)} (5.6) 



where 



a = b (z (z(l - v)K\ k) + i(l - v)^j , k 



a 
b' 



(5.7) 



Exploiting the explicit expression for G(p) found in [§] one can derive the following 
useful relation between G(—l) and G(—l) 



G(-l) 



Vksn(iuK',k)G(-l] 



ab) eG(-l). 



(5i 



Now, let us for a moment go back to our original model With a quadratic po- 

tential the model is invariant under the transformation M — > — M and the eigenvalues 
of the matrix M must hence live on an interval of the type [—a, a}. This means that 
for the support [a, b] of the eigenvalue distribution of the matrix X, defined in (|3.1|), 
we have b — -. Exploiting (|5.8|) and setting b = 1/a and we get from (|5.5|) and ( p.6|) 



d_ 

dp 



G{p) 



d 2 

dp- 2G(p) 



P =-i 



P =-i 



(1 - a 2 ) 2 



a^-e 2 + l | +(<\<~-<n)}G(-l). 



1 — a 



2\2 



1 - ae + 3ea + 2e 2 - a 2 - a^- } G(-l 



(5.9) 
(5.10) 
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Finally inserting (|5.9p and (|5.10|) into (|5.3|) and (|5.4|) both equations reduce to 



2d 1 

(2 - n) - -^—^ - {2 + 2ea + e~ 2 - a 2 } = 0. (5.11) 
5.2 The string susceptibility 

In this section we will determine the quantity ^T 3 ^ which we will make use of later 
when investigating the critical behaviour of the model. Here F stands for the genus 
zero contribution to the free energy of our model. The quantity ^T 3 ^ is related to 
the string susceptibility U(T) = ^ (T 2 F) by 

£_ i d _F\ _ dim 

dT 2 \ dT)' 1 dT ■ {bAI) 
By direct computation we find 

dT 2\N \1 + XJ I 2 Id 2iri \1+uj) K J K J 

Multiplying by T and differentiating once more gives 

d ,dF 1 f duj /1-uA 2 d /rTnTTr/ xx 

Now, it follows from ( |3.5| ) that (TW(p)) fulfills the following equation 

^ (TW(p + i0)) + A (TW(p - i0)) + n^(3W(-p)) = V» (5.15) 

where 

^(p) = (2-n)-i-. (5.16) 
1 +p 

Furthermore we obviously have for the asymptotic behaviour 

A(2W(p))~-, as p^oo (5.17) 

and as regards the analyticity structure, (XW(p)) must be analytic in the complex 
plane outside the support of the eigenvalue distribution and behave as 

4-(TW(v)) ~ (p-a)- l/2 , (p-b)- 1/2 for p -> a, b (5.18) 
aT 

Let us introduce the following notation 

W T (p) = -^{TW{p)) (5.19) 
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and let us split Wt(p) in a regular part, W^(p), and a singular part, W^{p), i.e. 

W T (p) = W£(p) - W£(p) (5.20) 
where W^(p) does not have any cut. Then we have from ( p. 15] ) 

2V\p) - nV'(-p) 



Wr(p) 



1 



n 



n 



P 



(5.21) 



4 — n 2 1 — p 2 

The singular part of Wt(p) is a solution of the homogeneous saddle point equation and 
as any other such solution has a parametrization of the form ( [4.3| ). Since WJ(p) has 
poles at p = ±1, V^Kp) must likewise have poles here because the full function Wt(p) 
should be analytic outside the support of the eigenvalue distribution. Therefore we can 
write 

W t(p) = Y^-i {A(p 2 )G(p)+pB(p 2 )G(p)} (5.22) 

where A(p 2 ) and B(p) 2 are now entire functions. From the requirement ( |5.17|) on the 
asymptotic behaviour and the expression (|5.21| ) for W^{p) one can conclude that A{p 2 ) 
and B(p 2 ) must be constants. Hence we have 



W T (p) 



n 



p + AG{p) + pBG{p) 



(5.23) 



p z 1 2 + n 

and the constants A and B are determined by the requirement that the poles at p = ±1 
should vanish, i.e. 

1 + AG{1) + BG{1) = 0, 



2 + n 
2 — n 

— + 1 + AG(-1)-BG(-1) 
2 + n 



0. 



The solution reads 



1 f 2nG(-l) -4(5(1) 1 
2 + n \G(1)G(-1) + G(-1)G(1)J ' 

1 [ 2wG(-l)+4G(l) 1 
2+^ +G(-1)G(1) J ' 

Going back to ( |5.14| ) we can write 



.4 



B 



(5.24) 
(5.25) 

(5.26) 
(5.27) 



rp6 

dT dT 



duo ( 1 — to 



2 Jox 2iri V 1 
dio 1 



LO 



LO 



LO ~ 



{AG(uj) + ujBG(uj)} 



2 /oi 2m (1+u) 3 

G( P ) 



{aG(uj) + ujBG(uj)} 



A 
~2 



dp 2 



d_ 

dp 



B 
~2 



' d d 2 \~ 
^ dp dp 2 J 



P =-i 



2 + n{l -a 2 ) 2 



[2 + e 2 - a 2 + 2ae} 



(5.2* 
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5.3 The critical behaviour 



As argued earlier our model becomes singular as a — > (cf. to sections ^ and |3|). Below 
we will investigate the nature of the critical behaviour associated with this singularity. 
In analogy with what was the case for the ordinary 0(n) model on a random lattice 
the present model only has a well defined scaling behaviour as a — > if n e [—2, 2] and 
we will restrict ourselves to considering this range of n values. One might also try to 
look for a critical point associated with a — > 1 (cf. to equation ( 5.11Q ), i.e. with the 



physical and the unphysical cut degenerating to the two points +1 and —1. Due to 
the analogy with the Penner potential |H| one might expect that having a = 1 (apart 



from at c = 0) is possible only for a particular value of n. (If the analogy were perfect 
it would be n = 1). However, we find that the equation a = 1 only has the trivial 
solution c = regardless of the value of n. 

5.3.1 The case n e] -2,2[ 

Let us consider the singular behaviour which occurs as a — > 0. First, let us fix n and 
determine the critical value of T as a function of n. By analysing the k — > limit of 
the various elliptic functions which enter the equation (|5.11| ) one concludes that in the 
limit a — > the dominant term in the curly bracket is a 2 and that 

1 

as a — > 0. (5.29) 

a 

Hence the critical value, T* of T is given by 

x 2z/ 2 . . 

(2 - n) - — = (5.30) 

or 

v 2 



^~ 2sin>7r/2)- (5 ' 31) 
In particular we see that T* is always positive and greater than 2/tt 2 . For n = 1 we get 
T* = |. In reference [|[ a numerical determination of this quantity gave ^- = 4.504. 



The next to leading order term in the curly bracket in equation (|5.11|) comes from 
the term 2ea which behaves as 

2ea ~ q {1 - u)/2 ~ ~ [^] " , a - 0. (5.32) 

(X Qj \ 4r / 

From this we conclude that 

T„ - T ~ a 2 " 2 ^. (5.33) 

Now, let us take a look at the expression ( |5.28| ) for ^T 3 ^. Here the leading order 
contribution comes from the term a 2 and is of order a . The next to leading order 
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term comes from 2ae and is of order a 2u (cf. to equation ( J5.32 )). Bearing in mind the 
relation ( |5.12| ) we get using ( [5.33D 



U(T) ~ (T* -T)t^. 



This means that 



5.3.2 The cases n = ±2 



1 - 1/ 



(5.34) 



(5.35) 



For n = ±2 the relations (|5.11|) and ( |5.28| ) contain divergent terms. However, the 
limits n — > ±2 of these relations are well defined. 



The case n = +2: Taking the limit n — > 2 in ( |5.11| ) one arrives at the following 
equation 

vr 2 - 75-^3- {(# + a 2 K'f -a 2 (l + a 2 ) V 2 } = 0. (5.36) 

This reproduces the result ( |5.31|) that T* = 2/7T 2 for n — 2. In the limit a — > the next 
to leading order contribution in the curly bracket comes from the term (aK') 2 which 
behaves as (a log a) 2 . This gives 



T* — T ~ (a log a) 
Furthermore, in the limit n — > 2 the relation ( |5.28| ) reads 



dT dT 4(1- a 2 )^ 



:i + a 2 )2 _ 4 z T7 



E' 



Letting a — > we get 



d »dF 1 

rp6 



(5.37) 



(5.38) 



(5.39) 



dT dT log a 

The results ( |5.37| ) and ( |5.39| ) coincide with those for the ordinary 0(2) model on a 
random lattice. 



The case n = —2: For n = —2 the relation (5.11) reduces to 

2 -(T^M (1+<i2)2 - 4 f} = 



(5.40) 



which in accordance with (|5.31|) gives that T* = i. Furthermore it follows that in the 
limit a — > 

r„-T~-^-. (5.41) 



log a 
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The relation (J5.28) takes the following form when n = —2 



In the limit a->0we find 

_lT 3 ^~-(aloga) 2 . (5.43) 

We note that the results (|5.41|) and (|5.43|) do not coincide with those of the ordinary 
0{— 2) model on a random lattice which (for gaussian potential) does not have any 
singular points. 

6 Conclusion and outlook 

We have solved exactly a hermitian (n + l)-matrix model with plaquette interaction. 
For n g] — 2, 2] the model was shown to belong to the same universality class as the 
0(n) model on a random lattice. In particular this result confirms the speculation of 
reference || that the critical point of the model ( |1.3|) describes the same physics as the 
critical point of the 0(1) model on a random lattice. Using equation ( 5.28Q it is easy to 



see that the plaquette model has no singular points (with T finite) for n < —2 and that 
for n > 2 the points given by vK' = 2mK, where v = iu, are singular. We expect that 
in analogy with the ordinary 0(n) model, the solution of the plaquette model breaks 
down at the first of these points, vK' = 2K, and that the critical index 7 str takes the 
value +| at this singularity. Although our model is much simpler than general lattice 
gauge models and matrix models generating Meander numbers our results may be taken 
as an indication that elliptic functions might provide a convenient parametrization of 
such models. 

Our solution of the plaquette model contains the solution of a certain three colour 



problem on a random lattice ||. The classical three colour problem due to Baxter [O 
consists in enumerating all possible ways of colouring with three different colours the 
links of a 2D regular three coordinated lattice in such a way that no two links which 
meet at the same vertex carry the same colour. The problem can also be understood as 



the problem of counting all possible foldings of the 2D regular triangulated lattice |12 |. 
Obviously the partition function (|1.3| ) (with V(A) = \A 2 ) generates random lattices 
with links of three different colours where no two links radiating from the same vertex 
have the same colour. Due to the matrix nature of the fields, however, in the present 
case the cyclic order of the three colours around a vertex will always be the same. 
In order to lift this constraint we would have to introduce two interaction vertices, 
tr ABC and tr ACB. The quartic interaction term in the resulting two matrix model 
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would then look like c (tr ABAB + tr A 2 B 2 ) |5|]. Unfortunately an exact solution of a 
model with this type of interaction is still lacking. Let us mention in this connection 
that a somewhat similar interaction term, namely c ( tr ABAB + 2 tr A 2 B 2 ) appears in 
a matrix model describing an Ising spin system living on the vertices of a randomly 



quadrangulated surface [13 
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Appendix A The Virasoro constraints 

Like the 0(n) model on a random lattice our model can be understood as a deformation 
of the one-matrix (n=0) model and obeys a set of Virasoro constraints obtainable from 
those of the n = model by a canonical transformation [ I4"| . Let us rewrite ( j2.1| ) as 

Z n = y rfMexpjtr ^t,M t |exp|-^-(trM fc ) j. (A.l) 

Introducing a differential operator, H, by 

-ssi(£)' 

we have 

Z n = e n Z . (A.3) 
The n = model obeys the Virasoro constraints L m Z = 0, m > — 1 where 

moo 00 o 

*m = £^r-+E***er- (A,4) 

k—0 ub k ub m—k k=0 ub m+k 

and 

J-Z = NZ . (A.5) 
The general model obeys the Virasoro constraints L m Z n = 0, m > — 1 with 

L m = e H L m e- H = L m + njr ^r^— ■ (A.6) 

fc=0 ot m +k 
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